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Since Schwarzshild discovered the point-mass solution to Einstein's equations that bears his name, 
many equivalent forms of the metric have been catalogued. Using an elementary coordinate trans- 
formation, we derive the most general form for the stationary, spherically-symmetric vacuum metric, 
which contains one free function. Different choices for the function correspond to common expres- 
sions for the line element. From the general metric, we obtain particle and photon trajectories, and 
use them to specify several time coordinates adapted to physical situations. The most general form 
of the metric is only slightly more complicated than the Schwarzschild form, which argues effectively 
for teaching the general line element in place of the diagonal metric. 
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I. INTRODUCTION 

For many decades, general relativity was considered 
a highly abstract subject, the purview of sophisticated 
mathematical physics far removed from the undergrad- 
uate curriculum. This situation has changed in recent 
years as general relativity has matured into a true ex- 
perimental science: gravitational lensing now routinely 
probes dark matter in distant galaxiesi, gravitational 
wave detection thresholds are steadily marching into the 
range of realistic sources^, and the Global Positioning 
System^ must compensate for gravitational time delays. 
In parallel to this evolving status of general relativity, a 
number of excellent textbooks at the undergraduate level 
have been published in recent years^Si^, demonstrating 
that teaching general relativity to physics majors is not 
only possible but exciting and relevant. 

One of the elementary results of general relativity is 
the spacetime arising from a point mass M, which is 
spherically symmetric and stationary, or independent of 
the time coordinate. Close to the mass, the spacetime 
represents a black hole with an event horizon, while at 
a large distance it reproduces the results of Newtonian 
gravity. This spacetime is most commonly represented 
by the Schwarzschild metric 



g^dx^dx 1 ' = x{ r )dt 2 



Xiry'dr 2 



2M 

X(r) = 1 , 

r 



r 2 dil 2 , (la) 



(lb) 



where we have used units c = G = 1 and abbreviated 



dn 2 



d0 2 



sin 



Historically, this was the first 



widely-used form of the metric^, and virtually all texts 
begin (and often end) the discussion of spherically sym- 
metric metrics with Eq. JQ). Formally, this form of the 
metric is in one sense the simplest possible, since it is the 
only diagonal form of the metric, and this property re- 
sults in simplified calculations (e.g., the smallest possible 
number of non-zero connection coefficients) . 

However, metrics in general relativity are not physi- 
cally measurable quantities, and numerous metrics em- 
ploying different sets of coordinates can be used to de- 
scribe the same physical spacetime. The usual form for 



the Schwarzschild metric possesses some peculiar prop- 
erties which make it far from ideal pedagogically. Most 
notably, it contains an unphysical coordinate singularity 
at the event horizon r — 2M, the interpretation of which 
is confusing enough that its true nature was debated for 
many years. A second more subtle point is that the met- 
ric is the same under time-reversal: a particle traveling 
forward in time is equivalent to one traveling backward in 
time, but physically particles can cross the event horizon 
in only one direction. 

Here we perform a general coordinate transformation 
of Eq^to obtain the most general form of the stationary, 
spherically-symmetric metric: 



ds 2 



X(r)du 2 + 2B(r)dudr 

- B(rf)dr 2 



r 2 dn 2 . 



(2) 



This form contains one free function of radius, B(r), 
which explicitly displays the coordinate freedom in the 
metric. Through various choices for B(r), all of the 
familiar (and not-so-familiar) forms of the metric can 
be obtained, including the Schwarzschild form and the 
Eddington-Finkelstein metric often used to investigate 
the physical nature of the event horizon. We advocate 
the use of Eq. @ as a much clearer and pedagogically 
versatile form of the metric for a number of reasons: 

1. Eq. J5J) emphasizes the fact that the same physical 
spacetime is described by many different metrics; 

2. It unifies the disparate expressions for the metric in 
the general relativity literature into a single simple 
form; 

3. It naturally suggests variables describing motion of 
particles which are closely tied to physically mea- 
surable quantities, and investigating the reduction 
to different metric forms clarifies the spurious na- 
ture of the singularity in Eq. and 

4. Mathematical manipulations involving Eq. (0) are 
only slightly more complex than those using the 
Schwarzschild form. 
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While Eq. (JIJ may be the simplest form of the metric, it 
does not give the simplest or most physically transparent 
forms for equations of particle motions in the spacetime. 

The rest of this paper derives the general form of 
the metric, demonstrates its reduction to various specific 
forms, and discusses the corresponding expressions for 
the geodesic paths followed by light and particles mov- 
ing in the radial direction. In the next Section, Eq. © is 
derived from the familiar Eq. via a straightforward co- 
ordinate transformation. Section IIIII establishes expres- 
sions for radial geodesies in the general case, and Sec. II VI 
examines how the geodesic expressions are simplified for 
specific choices of the arbitrary function B(r), includ- 
ing all of the most commonly used forms of the metric. 
Three brief appendices cover relevant but somewhat more 
technical material: Appendix lAl gives the connection and 
Ricci curvature components associated with Eq. © for 
convenience in doing calculations; Appendix iBl discusses 
why a change in time coordinate does not affect the grav- 
itational redshift; and Appendix IU1 expresses the general 
line element in vierbein form. The derivations and argu- 
ments in the body of the paper require no mathematics 
beyond elementary calculus, and as such are suitable for 
inclusion in any undergraduate general relativity class. 



II. THE GENERAL LINE ELEMENT 

In the Schwarzschild form of the metric, Eq. QJ, 
the time direction is orthogonal to all other directions 
(that is, gok = when k ^ 0) so the metric is called 
static; Eq. is the unique static spherically-symmetric 
vacuum line-element. The radial coordinate r in the 
Schwarzschild geometry has a clear physical interpreta- 
tion: it is the proper distance measured by a massive 
particle falling into a black holei The time coordinate t 
is more problematic. The proper time for a particle (ei- 
ther massive or massless) falling through the event hori- 
zon is finite, but the coordinate time t diverges — it takes 
an infinite amount of time in the Schwarzschild system 
of coordinates to fall into a black holei Thus the singu- 
larity of Eq. at the point r = 2M is not physical, and 
can be transformed away with a different choice of coor- 
dinates. We are interested in forms of the metric that 
are stationary but are not necessarily diagonal. 

A transformation of the time coordinate of the type 

u = at + P(r), (3) 

for some function fj and constant a, ensures that the 
metric remains stationary, while keeping the spherical 
symmetry manifest. In other words, if a and remain 
unspecified, we obtain the most general stationary vac- 
uum line element in spherical coordinates. Without loss 
of generality we take a = 1; a simple rescaling of the 
time variable is equivalent to choosing a different time 



unit. Since u is a function of t and r, its differential is 

, du du dp 
du = —at + -^—dr = at + — dr 
ot or or 

du 2 = dt 2 + 2P'(r)dtdr + P'{rfdr 2 . (4) 

Rather than going through the complicated algebra of 
solving for dt and substituting into Eq. (JJJ, note that a 
general stationary line element must have the form 

ds 2 = A(r)du 2 + 2B(r)dudr - C{r)dr 2 - r 2 dfl 2 ; 

setting this equal to Eq. and using the relations in 
Eq. gj, 

A(r) = x(r), (5a) 
B(r) = - X (r)P'(r), (5b) 
C(r) = X (r)P'(r) 2 + 2B(r)p'(r)+x(r)' 1 . (5c) 

We can use Eq. (|FSb|) to eliminate p' = —B(r)/x(f), which 
leads to a condition on B and C via Eq. (|5*c|) : 

C(r) X (r) + B(r) 2 = 1, (6) 

and Eq. (J2J) follows. As Appendix lAlshows explicitly, the 
Einstein field equations do not specify B(r), so we will 
call Eq. J3J) the general line element for the spherically- 
symmetric vacuum equations. We are free to choose 
B{r): for example, the Schwarzschild form of the met- 
ric is obtained from the general line element by letting 
B(r) — 0. We will discuss other coordinate choices below. 
Appendix IBl examines how a change in time coordinate, 
specified by B(r) and x( r )> leaves some of the standard 
tests of general relativity unaffected. 

Of course, it is also possible to change the spatial co- 
ordinates (r, 8, <ft), but it is nice to leave them as written 
for two reasons: spherical symmetry is explicit, and as 
mentioned in the Introduction, the radial coordinate r is 
a measure of the proper distance measured by a massive 
particle falling radially into the black hole. 

III. RADIAL GEODESICS 

Since free particles (including photons) follow 
geodesies in general relativity, the Lagrangian for free 
motion can be written simply as the squared norm of the 
velocity x = dx/dX with respect to some parameter Ai 

2C(x, x, A) = g llv (x)x ,J, x u = xu 2 + 2Biir 

- x _1 (l - B 2 )f 2 - r 2 {9 2 + sin 2 d<j) 2 ). (7) 

The tests of general relativity that depend on the angular 
variables (light deflection, for example) do not explicitly 
involve the time coordinate, so we will look strictly at 
radial geodesies by letting 8 — ir/2 and <f> = 0. 

The magnitude of a geodesic vector (and hence the 
Lagrangian itself) is a constant of motion, so we have 

g liV x ih x v = xu 2 + 2Biir - x _1 (l - B 2 )f 2 = k (8) 
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where k = 1 for massive particles and k = for photons. 
The time coordinate u is cyclic, so we know that 



— — \u + Br — e 
au 



(9) 



is constant, with e related to the energy of the particle 
(with respect to an observer at infinity) as follows: 

{E/m massive particles 

E/Vh photons 

Solving Eq. © for u and substituting it into Eq. © 
allows us to solve for f: 



these examples are familiar, others less so. Our central 
conceptual point is that the general form of the met- 
ric, Eq. (J2J), allows us to choose a specific system of co- 
ordinates which leads to simple and physically intuitive 
forms for the geodesies. Since the geodesies, unlike the 
metric, represent physically observable quantities, pick- 
ing a coordinate system based on the properties of the 
geodesies instead of the metric makes sense, as long as 
the coordinate system does not lead to unduly complex 
computations. 



A. Two Simple Systems 



r = e — KX- 

Since this is independent of the choice of B, the radial 
geodesic equation is as well: 



r + nM/r 2 



0. 



(10) 



which is the Newtonian expression (upon restoration of 
G) — although r is not the flat spacetime radial coordi- 
nate, and when k = 1, A is not the coordinate time but 
the proper time shown by clocks moving with the massive 
particle along the geodesioi. 

Using our expression for r in Eq. @ allows us to write 
the velocities with respect to the proper time: 



w\J e 2 — kx 

{e - wBy/e 2 - kx} 



(11a) 
(lib) 



where w = +1 for outgoing geodesies, while w = — 1 for 
incoming. From these, we can find the coordinate time 
as a function of radius, 



u(r) 







[±dr = 




I r j 





w 



X I y/1 - «x/e 2 



B 



(12) 



which tells us whether a particle crosses the event hori- 
zon in finite coordinate time, for a particular choice of 
B(r). In the Schwarzschild case B = 0, for example, 
the time required diverges logarithmically for both val- 
ues of ui Once B is chosen, two parameters dictate the 
type of geodesic: k = or 1, which specifies lightlikc 
(photons) versus timelike (massive particles) geodesies, 
and e, which determines the energy of the particle. In 
the case of null geodesies (k = 0), the photon energy is 
irrelevant: 



u(r) = f — (w-B). 
J X 



(13) 



IV. 



COORDINATE SYSTEMS FROM 
GEODESICS 



Equation (|13|) gives the coordinate time as a function of 
radius for a null geodesic. Future-pointing paths of light 
in special relativity are straight lines defined by it(r) = 
—r + const.; to replicate this, we let B = x~ 1 = —2M/r, 
which means that 



u E F{r) = 



-r + u a 
r + AM lnlr 



2MI + u 



incoming 
outgoing 



where uq is a constant of integration. This is one ver- 
sion of the familiar Eddington-Finkelstein coordinates; 
ingoing geodesies are straight lines, while outgoing must 
"climb out" of the gravity well formed by the black hole. 
Nothing starting within the event horizon can get out, 
since the "outward" geodesies starting from r < 2M ac- 
tually fall inward as well. 

The line element for this choice of coordinates is 



ds 2 EF — du 2 



dr z - r^dQ, 
(77^ - 2Vl IM Q dx»dx v 



— - (du — dr) 2 
r 



where in the second line we have written the metric as 
the sum of the flat metric (in spherical coordinates) and 
a contribution from a null vector 

H = (1,-1,0,0). 

Metrics of this type are often called Kerr-Schild metrics, 
and in their more general form are used to find the ro- 
tating black hole solution (known as the Kerr solution; 
see Ref. 01 f° r more details). 

A related class of coordinates can be obtained by set- 
ting the proper time (for a massive particle) or affine 
parameter (for a photon) equal to the coordinate time 
for an inbound geodesic. This means we must specify in 
advance the energy e' of the particle, and whether it is 
massless or massive by setting the parameter k': 



u = 1 



X 



[l + B^l-n'x/e 12 ] 



Selecting the retarded solution, we find 



In the following subsections, we will examine choices 
of B(r) that simplify the geodesic expressions. Some of 



B = 



1 



£-1 



(14) 
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where 



P ; 



k /e 



is the metric parameter that determines the time co- 
ordinate associated with a particular energy. When 
k' = p = 0, B = x/e' - 1; if e' -> 1 (so that |r| = 1), 
we get back the Eddington-Finkelstein time coordinate. 
Thus, despite the appearance of both p and k' in Eq. 1)140. 
they are noi independent parameters. Since we have al- 
ready discussed the Eddington-Finkelstein case, we will 
consider only k' = 1 and p ^ for the remainder of the 
section. 

The line element for this geometry with k' = 1 is 



ds 2 G T = Xdu 2 



dudr 



V 1 -px 

2^p-p{l + X ) jj> 
1 - PX 



-dr 2 - r 2 dft 2 (15) 



where the subscript "GT" indicates geodesic time. For 
an arbitrary geodesic, the trajectory is described by 



M r ) = I ~ 



xVp- i 



X { y/l - Kx/e 2 V 1 - PX J 

Plots for different values of p and e are shown in Fig.^i. 
When e = e' (and k = 1), the time for an infalling 
geodesic takes the simple form 



u p (r) = dr 



P 



i-px 



which reduces to the Eddington-Finkelstein case in the 
limit p — > 0. 

The relative complexity of the line element in Eq. i|15|) 
may be why it is not more widely known — except in the 
specific case p = 1, which is the Painleve-Gullstrand met- 
ric (see Ref. \lHt) . independently discovered by Lasenby, 
Doran, and GulU£: 



dsp G = xdu 



2M 



-dudr — dr 



r 2 dn 2 . 



(16) 



We will see this metric again as a special case in the 
following section. 



B. Relationship Between u and Proper Time 

In the previous section, we have chosen coordinate 
times u which lead to physically intuitive forms for ra- 
dial geodesies, specified by u(r). In one case we explic- 
itly equate the coordinate time with the proper time for a 
massive particle moving along a radial geodesic, Eq. (|15fl . 
More generally, we can derive the proper time u for any 
choice of coordinate time u (or equivalently, the function 
B{r)). Starting from Eq. (JSJ), we multiply both sides by 
dX/dr and rearrange somewhat: 

4-(«sH5-Hff (17) 



Applying Eq. (|llbjl to the first term of the right-hand 
side eliminates du/dX 



du \ u du du 

X- r +2B) - = — - — w B^\-px 

QiT / c ' ! ' ' 



dr dr 
IB 

+ (1-wByJl-px) 

X 



du 
dr 



B(2 X - 1 -1) {l-wB^rTpx, 

where p = k' /e' 2 as before, and we have used the identity 
du/dr = u/f. Substituting Eq. (|lla[> and integrating 
both sides of Eq. I|17|) over r yields 



e'pX 



dr 



B 



V 1 - px) 



(18) 



which is a coordinate transform as in Eq. For massive 
particles with p ^ 0, Eq. I|18|) gives the proper time A in 
terms of the coordinate time u; when p = 0, no proper 
time can be defined but the integral is the same as the 
one in Eq. (JT3J). 

For special choices of B, the integral in Eq. (|18|) takes 
a simpler form. The obvious example is 



B = -^/T~p^, (19) 
which gives rise to the Lake-Martel-PoissoniiiH metric 



ds 2 — xdu 2 — 2\J\ — pxdudr — pdr 2 



2 dn 2 



Different values of p yield several well-studied solutions: 
when p = 0, we obtain the more familiar form of the re- 
tarded Eddington-Finkelstein metrioi (in which infalling 
null geodesies are lines of constant u) , while p = 1 yields 
the Painleve-Gullstrand metric of Eq. ((ToT) . We may also 
treat p in general as a positive parameter (which can be 
associated with the energy of a particle at infinit}iii)in 
which case r can range from zero to infinity for < p < 1 , 
and 

r<2Mp/(p-l) for p > 1. 

The latter case is proportional to the Gautreau- Hoffmann 
coordinate system^, with the initial radial value defined 
by Ri = 2Mp/(p—l) in our notation. Note that Eq. fl^l 
is equivalent to Eq. (|14fl when p = e' = 1; in this limit, 
by construction, the proper time is equal to the coordi- 
nate time. This case is equivalent to setting u — e'p and 
repeating the analysis of Section liV Al Plots for different 
values of p and e are given in Fig. \I]p. 



V. DISCUSSION 

The general line element provides a clear and straight- 
forward way to pick specific time coordinates for the 
spherically symmetric, stationary spacetime. Eq. (J2J) is 
a simpler starting point than the diagonal Schwarzschild 







2M 





r 



2M 



FIG. 1: a. Plot of the geodesic time parameter ugt(t) for various values of p. The legend indicates different values for the 
energy parameter e, as well as the null geodesic, b. Plot of the Lake-Martel-Poisson time parameter ulmp{t) under the same 
conditions as the geodesic time. Note that the times are identical (ugt = ulmp) for p = 1. 



metric for determining the two forms of the Eddington- 
Finkelstein coordinate system, and provides both the 
geodesic time parameter and the Lake-Martel-Poisson 
coordinates directly, as opposed to performing a rather 
complicated coordinate transformation. In fact, the 
geodesic time line element is new to us, and obtaining 
it directly from the Schwarzschild solution involves a dif- 
ficult integral. 

We have emphasized stationary forms of the metric; 
however, it is possible to "maximally continue" any of 
these coordinate systems to make them diagonal. This is 
done by using 

u = t + (3{r) and v = t — 0(r) 

in place of t and r, and rewriting the line element in 
the new coordinates. The familiar example of this is the 
Kruskal-Szekeres extension, which uses the Eddington- 
Finkelstein coordinates as its starting point. However, 
all maximally-extended solutions include a non-physical 
region of spacetime (the "white hole" problem), and this 
poses as serious an interpretational difficulty as the co- 
ordinate singularity problem at the event horizon. 

From a pedagogical point of view, we emphasize that 
Eq. (J2J contains virtually all other discussed coordinate 
systems as special cases. The general metric is only 
slightly more complicated to use in calculations than 
the diagonal Schwarzschild case; for convenience we list 
the connection and Ricci tensor components explicitly 



in Appendix 1X1 For reasons of pedagogical and concep- 
tual clarity, We advocate presenting Eq. as the line 
element describing the spherically symmetric and sta- 
tionary vacuum spacetime when teaching general rela- 
tivity, especially at the undergraduate level; the familiar 
Schwarzschild form of the metric can be easily obtained 
as a special case. Appendix lUlshows how to derive Eq. 
in a straightforward way which is no more complicated 
than usual derivations of the Schwarzschild metric. 
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APPENDIX A: SOLUTION OF EINSTEIN'S 
EQUATIONS 

The Christoffel connection coefficients are found using 
the standard formula 

r ^ = \$ Xk (d^g KU + d v g Kll - d K g^) . 



6 



If we take Eq. J5J) to be an ansatz for solving Einstein's 
equations, with unknown functions x( r ) an d B(r), the 
nonzero connection coefficients are 



1 oo — 

1 oi — 

1 ii — 

1 22 — 

1 oo — 

1 22 — 

r 2 - 

1 12 — 

r 2 - 

1 33 — 

r 3 - 

1 23 — 



-Y l 0l = B X '/2 
-rii = C X '/2 
CB' - BC'/2 
T° 33 esc 2 6 = Br 
XX72 



r 3 - r- 1 

1 13 — ' 

— sin cos 
cote 



where we have used the shorthand (see Eq. iJBJ) 

C = X - 1 (l-B 2 ) 

and x' — dx/ dr. These expressions are not complicated, 
but there are more of them than in the diagonal case. 
The Ricci tensor components, 

p n per _i_ -pp per pp per 

J^HU — U V L ^ a U a L fj,v ' 1 pa L pv 1 pv L pa' 

turn out to be quite simple: 



-Roo 
-R22 



X$/2r Rot = B<P/2r R n = -C$/2r 
i? 33 esc 2 6 = 



where 



$ = X "r + 2 X '. 



Thus, Einstein's vacuum equations = do not in- 
volve B(r), and are easily solved to give Eq. Ijlbfl . with 
M as an integration constant which is identified with the 
total mass by ensuring the Newtonian limit. 



For gravitational redshift, let = (1,0,0,0) be a 
timelike (Killing) vector in whatever coordinate system 
we desire. The "direction of time" at a given value of r 
is then indicated by 



(Bl) 



If we locate two stationary observers 0\ and O2 at R\ and 
R2, respectively, and send a pulse of light between them, 
we obtain a shift in frequency due to the presence of X - 
The path of the light is indicated by the null momentum 
vector k, whose time component gives the frequency: 



r=R 2 



The projection of the timelike vector £ along the geodesic 
is constant (although proof of this conjecture requires 
a discussion of Killing vector o 8 ! 14 , which is beyond the 
scope of this paper). Using Eq. i|Bl() . we obtain 



wi = k^e/Vxin) 

and a similar expression for u>2- Taking the ratio elimi- 
nates the dependence on k, leaving us with 



ui2 



' x{r2) 



(B2) 



The redshift formula is directly obtained from this ratio; 
since it doesn't depend on the choice of time coordinate, 
it holds for all variations on the general line element. 
Thus, the redshift is infinite if the emitter is located at 
n = 2M, even if the coordinate system allows for a cross- 
ing of the event horizon in finite time. 



APPENDIX C: VIERBEIN FORM OF THE 
GENERAL LINE ELEMENT 



APPENDIX B: NEW TIME VARIABLES AND 
THE "STANDARD TESTS" 

Most derivations of the measureable effects due to a 
planet, star, or black hole start from Eq. (JIJ; neverthe- 
less, changing time variables does not change the predic- 
tions of general relativity. We will not bother to rederive 
all the familiar expressions, but instead will emphasize 
their independence of particular coordinate assumptions. 
In fact, most tests of GR (bending of light, perihelion ad- 
vance of an object in orbit, etc.) depend on the radial 
velocity with respect to proper time, which is a general- 
ization of Eq. (|llal) . Since this equation does not involve 
u in any way, we don't need to discuss these tests here. 
The exceptions are time delay and gravitational redshift; 
demonstrating that the former is independent of time co- 
ordinate choice requires coverage of non-radial motion, so 
we will not discuss it here. 



The general line element expresses all the metric de- 
grees of freedom that are possible by changing the time 
coordinate. However, some degrees of freedom lie "be- 
neath" the metric; these are often expressed as objects 
called vierbeins h®, which are like a "square-root" of the 
metric. The equivalence principle allows us to assign an 
orthonormal frame of vectors (whose metric is just the 
flat spacetime metric r) a b) to each point of spacetime, and 
the vierbeins take us between the orthonormal frame and 
the full metric of GR: 

9fiv = h®hlr) a b -> \J-dctg = det h. 

Writing things out in terms of the vierbeins gives us not 
only the general line element, but also lets us express 
the isometries — coordinate transformations that do not 
change the form of the metric. 

We write in matrix form the vierbeins associated with 
the most general spherically symmetric and stationary 
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spacetime as 



and 



so that 

ds 2 = (g 2 



( 9i(r) /a(r) 

92 (r) h(r) 

r 

V rsinfl 



g\)du 2 + 2(/i#2 - hgi)dudr 



(f?-fl)dr 2 -r 2 d{l 2 . (CI) 



(This is the way we first obtained the general line 
element.) Eq. J2J follows from the identifications 



X = 9\ 



92, 



B = fw2-f 2 9i, C = fl-f%, 



fm - 1292 = l. 



(C2) 



Einstein's equations then give Eq. (|lb|) . and Eq. (|C2|) 
leads directly to Eq. J(|J|, so at first glance it appears that 
these two conditions on the four free functions leave two 
free functions to be specified. Nevertheless, the condi- 
tion C\ + B 2 = 1 arises from f-ygi — f%gi = 1, so that 
the metric still only contains one free function. The two 
free vierbein functions are both required to determine the 
metric, so individually they represent isometric degrees 
of freedom — a given expression for B corresponds to a 
great many choices of vierbein fields. 
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